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Abstract 

Using unitary transformations, we express the Kondo lattice Hamiltonian in terms of 
fermionic operators that annihilate the ground state of the interacting system and that rep- 
resent the best possible approximations to the actual charged excitations. In this way, we 
obtain an effective Hamiltonian which, for small couplings, consists in a kinetic term for 
conduction electrons and holes, an RKKY-like term, and a renormalized Kondo interaction. 
The physical picture of the system implied by this formalism is that of a vacuum state 
consisting in a background of RKKY-induced spin correlations, where two kinds of elemen- 
tary modes can be excited: Soft neutral modes associated with deformations of the spin 
liquid, which lead to very large low-temperature values of the heat capacity and magnetic 
susceptibility, and charged modes corresponding to the excitation of electrons and holes in 
the system. Using the translational and spin rotational symmetries, we construct a simple 
ansatz to determine the charged excitations neglecting the effects of the spin correlations. 
Apart from the 'normal', uncorrelated states, we find strongly correlated charged modes 
involving soft electrons (or holes) and spin fluctuations, which strongly renormalize the low- 
energy charged spectrum, and whose energy becomes negative beyond a critical coupling, 
signaling a vacuum instability and a transition to a new phase. 



PACS: 71.27.+a, 75.30.Mb, 75.20.Hr. 
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I. INTRODUCTION 



Despite many years of intense experimental and theoretical studies of heavy-fermion 
(HF) systemsjlhi their physical picture still remains controversial. It is generally believed 
that both HF metals and the so-called mixed-valence compounds are well-represented by 
the Anderson lattice model (ALM) and that their phenomenologies strongly depend on the 
location of the f band with respect to the Fermi level. Ce-based HF compounds are usually 
regarded as an extreme case in which the deep f level and high on-site Coulomb repulsion 
strongly inhibit f-charge fluctuations. This has led many authors to regard the Kondo 
lattice model (KLM) as a canonical model for these compounds since it has been proved 
to be the limit of the ALM for vanishing f-charge fluctuations,! and it exclusively involves 
the degrees of freedom that are expected to be relevant to the low-energy phenomenology 
of these materials. 

Among the most important features characterizing this phenomenology are the very large 
values of the low-temperature magnetic susceptibility and linear coefficient of the specific 
heat, the appearance, in some cases, of antiferromagnetic order with strongly quenched local 
moments, and the transition to a superconducting phase in systems like CeCu2Si2- 

HF metals are often regarded as 'dense Kondo systems' and its properties qualitatively 
discussed extending to the lattice the results obtained for the singleii and two-impurityi 
cases. In the first of these cases (Kondo problem), a singlet state is formed at zero tempera- 
ture between the spins of the impurity and of the conduction electrons.! Since the low-energy 
singlet-triplet excitation associated with this (Abrikosov-Shul) resonance gives a contribu- 
tion to the specific heat,0 the formation of one of such resonances at each site of a Kondo 
lattice would account for the very high heat capacities measured in these systems. In the 
two-impurity case, however, we have an additional ingredient, which is the appearance of an 
induced, inter-site (RKKY) interaction. Thus, magnetic order in the Kondo lattice is often 
interpreted as due to a competition between the RKKY interaction (which tends to produce 
magnetic order) and the Kondo effect (which tends to quench the local moments) .0 
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These interpretations, however, have been criticized by Nozieres,0 who pointed out that 
only a very narrow layer of conduction electrons around the Fermi level are energetically 
capable of forming Kondo singlets, which implies that we can have only one 'efficient' electron 
per many local spins. Moreover, Nozieres argument implies that, if strongly correlated (SC) 
states between conduction electrons and local moments are formed in the system, they 
should not correspond to local Kondo singlets but to extended, collective states in which 
one electron couples to many local moments, namely, to a spin wave. The calculation of 
such SC states or 'magnetic polarons' is, in fact, one of the main objectives of this article. 

Another consequence of the Nozieres argument is that only the very-low-energy part of 
the electronic dispersion relation can be strongly renormalized and, obviously, this mass 
renormalization cannot account for a measured entropy at the characteristic Kondo lattice 
temperature (Tjcl) of the order of iVlnL, where L and N are the f-level degeneracy and 
number of lattice sites respectively. Since conduction electrons are the only source of electric 
charge in the KLM, charged fermions cannot be responsible for the huge masses implied by 
the low-temperature specific heat and magnetic susceptibility. 

In contrast to the slave-boson00 and variational^ approximations, which defend the 
existence of charged heavy particles, the above arguments have led Kagan, Kikoin, and 
Prokof'ev0 to suggest that the main contribution to the low-temperature thermodynamics 
of these materials comes from neutral Fermi excitations of spin origin. They showEl that this 
picture is in good agreement with the available experimental data while a charged-fermion- 
based approach seems to be in conflict with De Haas- Van Alphen experiments.^ 

In this paper, we introduce a mathematical formalism to study the KLM that leads to 
a physical picture of HF systems (summarized in the conlusions) which does not enter in 
conflict with the Nozieres exhaustion problem and which, among other things, supports the 
views presented in Ref. 16. 

The basic idea of this formalism is presented in Section II, and consists in carrying 
out a unitary transformation to determine a new set of canonical fermionic operators that 
annihilate the ground state of the interacting system. We find that, when the Kondo lattice 



Hamiltonian (KLH) is expressed in terms of these operators, it essentially consists, for small 
couplings, in a kinetic term for electrons and holes containing a gap in its dispersion relation, 
an RKKY term, and a Kondo interaction in which the terms involving an electron and a 
hole operator have been cancelled. 

Using this expression, we discuss magnetic order in Section III. We argue that, since all 
the terms of H that contain fermions annihilate the vacuum state determined in Section 
II, the RKKY interaction is the only one responsible for the dynamics of the spin lattice 
unless some electronic instability (which we indeed find for stronger couplings) takes place. 
Thus, we believe that the reasons for an eventual formation of a spin liquid should be 
seeked in the enhancement of spin fluctuations due to the strong frustration produced by 
the long-range RKKY interaction. We also argue, in agreement with Ref. 16, that the very 
large specific heat measured in these systems should be attributed to the enormous entropy 
increase associated to the thermal breakdown of soft spin correlations. 

We mentioned that the dispersion relation corresponding to the kinetic term for electrons 
and holes of the transformed Hamiltonian contains a gap. However, if no symmetry is broken 
and no bound states are formed, there is no physical reason for the appearance of a gap in 
the middle of the conduction band. In Section IV, we calculate the actual (uncorrelated) 
charged modes neglecting the effect of the spin correlations and show, to leading order, that 
no gap is actually opened at the Fermi level. 

As we mentioned when we discussed the Nozieres exhaustion problem, if SC states are 
formed in the system, they should be extended, collective states in which a conduction elec- 
tron or hole couples to a spin wave. In Section V, we look for such states. We argue that, 
to include the effect of the high-energy tail of the Kondo interaction, these modes should 
be constructed from 'optimal' fermionic operators that represent the best possible approx- 
imations to the actual (uncorrelated) charged modes. These new operators are determined 
with the aid of a second unitary transformation that leaves the ground state invariant and 
which is determined variationally. We find that the expression of H in terms of this new set 
of canonical operators essentially consists in a kinetic term containing a much smaller gap, 
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an RKKY term without any modification, and a residual Kondo interaction which is only 
effective at low energies. With these optimal operators and neglecting the spin correlations, 
we construct SC charged modes using an ansatz that exploits all the symmetries of the 
Hamiltonian. We find that SC modes of spin 3/2 and 1/2 can form, respectively, in ferro 
and antiferromagnetic Kondo lattices and that its formation is much more favored in the 
last case. We also find that for suffiently strong couplings, these modes must strongly renor- 
malize the low-energy charged spectrum and that, beyond a critical coupling, they would 
condensate in the ground state producing an instability and, consequently, a transition to 
a new phase. Since the coupling of spin waves to conduction electrons to form SC modes 
influences the spin correlations (which are the main source of low-energy entropy) , this phase 
transition should be accompanied by important changes in the specific heat. 

II. FIRST TRANSFORMATION: INTERACTING VACUUM STATE 

We shall study a KLM consisting in a conduction band coupled to a lattice of s = 1/2 
local moments by an exchange interaction. The Hamiltonian is given by 

H = 5Z £ kcL c ka + jJ2 S eiSfi, (1) 

k,a i 

where S e j and Sfi are the spin of the conduction electrons and the local spin at site i re- 
spectively. S e i = l/2Y^a,f3 c la a ai3Cif3, where c\ a are electronic operators in the site (Wannier) 
representation. If Rj is the position of the i site, these operators are given, for a lattice with 
N sites, by c\ a = N~ l l 2 S k e- !R ' k c[ Q . 

We will consider the case of a half-filled conduction band with a symmetric density 
of states. Although the symmetry between electrons and holes will be very convenient 
to simplify the formalism, the scenario presented in this paper is expected to hold in all 
situations around half filling. To simplify the notation, we shall take energy units in which 
the band width is normalized (W = 1). The final formulas that we will obtain can, of course, 
be written in arbitrary units by suitably reintroducing W using dimensional analysis. In ([!]), 
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we have already subtracted the chemical potential in and, thus, W = 1 and the electron- 
hole symmetry, imply that — 1/2 < < 1/2. Finally, to enable an analytic calculation of 
many quantities throughout this paper, we will assume a constant density of states D. D 
normalized and constant, and W = 1 imply D — 1. 

The fundamental operators Ck a , S/j are characterized by the following algebraic proper- 
ties: 



{cka,4'/3} = <5kk"W , {Cko,^} = 0, (2) 

: X^/^/i ' ( 3 ) 

7 

Sj, = S /4 , (4) 

S? 4 = 3/4, (5) 

[ck«,S/i] =0. (6) 

As it is well-known, Eqs. (|]) and (|J) can be satisfied by writing S/j = 1/2 Y, a ,p fia^apfifi, 
where f} a and fi a are creation and annihilation operators of / electrons satisfying canonical 
anticommutation relations. It can be proved (and is physically obvious) that the condition 
(§) is equivalent to demanding that each site be occupied by one and only one / electron. 
As it will soon become clear, instead of using the S /j operators, it is much more convenient 
to work with the following equivalent set of operators: 

s 1 , i = V2(S% + iS y fi )^V2fy il , 

so,i = 2Sfi = / iT /»t - fijii , (7) 
s^ i = V2(s%-iS%)^V2flf iV 

It is straightforward to see that these operators, acting on states containing a single / 
electron per lattice site, satisfy the following multiplication table: 



s 0,i s l,i — s l,i 


s Q,i s -l,i 


s l,i s 0,i = — s l,i 


s -l,i s Q,i 


Si :i S-i :i = 1 + S ,i 


S 2 
1)* 


S-tjSij = 1 - S ,i 


<? 2 
b a,i 



(8) 

Li " ''' 



£i , = 



1 



as well as the properties 

s I,t = S -l,i 5 S 0,i = S 0,i j (9) 

[s M , S ^]=0 Vz^j. (10) 



Eqs. (H)- (|10D can also be regarded as the defining conditions of the algebra of local 
moments since, as it can be readily checked, they imply Eqs. ©-(H)- 

The multiplication table (§) will be extremely useful since it will enable us to reduce 
products of si operators corresponding to the same lattice site. For example, 

S 0,5 S -1,7 S 1,5 S 1,7 = s l,5 — S 1,5 S 0,7 ■ (H) 

In order to exploit the symmetry of the Hamiltonian (JJ) under lattice translations, it is 
convenient to work with all operators in the momentum representation: 

<L = N- 1/2 E e^ k 4 , s l>k = IT 1 " E , (12) 

i i 

where k is always in the first Brillouin zone. The s^k operators satisfy the following equa- 
tions, which will later be of much use: 

[so,k, s±i,k'] = ±2A^~ 1/2 s±i ik+k / , 

[si,k, s-i.k'] = 2A^ 1/2 s ,k+k' , (13) 
[si,k, si tk >] = , 

4k = s -*,-k- (14) 

Another tool which we are going to use constantly is the development of a product of 
two Si operators in the momentum representation; here is an example: 

s-i, k si,k' = *(k + k') - N-^so^y + N- 1 E e^+^s^s^ . (15) 

Finally, in the momentum representation, the KLH (HI) is given by 
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J 



4AT1/2 



( c kt c k'T _ 4| c k'i) s o,k'-k + V^4.T c k'| s -i,k'-k + \/2cj ci c k ' T Si >k /_ k 



(16) 



where, as in many equations throughout this paper, sum over repeated indexes in each term 
of the equation is implicitly understood. 

We want to study this Hamiltonian for small values of J/W, which is the region char- 
acteristic of HF systems. In the absence of interaction, all the spin configurations will be 
degenerate, and the ground (or vacuum) state of the system |$ ) will consist in a completely 
disordered spin lattice plus a Fermi sea containing all electrons with negative energy. The 
charged excitations of this system are electrons and holes consisting, respectively, in adding 
or removing an electron from this Fermi sea. To simplify the notation, q (p) will always 
denote a wave vector located below (above) the Fermi surface (e q < 0, e p > 0), while k will 
be a wave vector without such restrictions. As it is obvious and well-known, the creation 
and annihilation operators corresponding to electrons and holes are given by 



(17) 



and the non-interacting vacuum state can be written as 

l$o>= E C^ a jl ai ...f NaN \%), (18) 

where |^o) is a state annihilated by e pa and hq a , and the coefficients C^y___ are given 
by C}$_ a = 2~ N l 2 e %RoL ^--- OL N , where R ai ...a N are random phases, as it corresponds to a 
completely disordered spin lattice. 

The KLH fll6|) is given, in terms of the operators (|I7D , by 



H — (J2q2 £ q) + £ p e Pa e pa £ q ^ qa ^ q a 
J 



4AT1/2 

J 



' e P T e p'T - e pl e P'l) s o,p'-p 



/_p + v / 2e pT e p ' i s_i i p/„ p + v/2e^ep/fSi )P /_p 



4AT1/2 
J 



(^qT^q'T ~ /l q|Vl) S 0,q'-q ~ V^h^k^ jS _l, q '- q - V^/l^ /l q /| S^q/ q 



4AT1/2 
J 



(eLhl - eLhL)s _q_ p + v / 2e_ T /iLs_ 1 _ q _ p + v^eL /i_ 



; p|"q| Sl -q-p 



4AT1/2 



(^qi e pT - ^qt e p|) s o,q+p + v2/i q -rep-[-si iq+ p + v2/i q ie p ;s_i )q+p 
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(19) 
(20) 
(21) 
(22) 
(23) 



This expression makes clear that |$o) is not an eigenstate (and, hence, not the ground 
state) of H when J ^ 0, because the terms (p2|) do not annihilate |$ ) but create additional 
electron-hole excitations accompanied by a spin fluctuation. 

To determine the elementary excitations of the interacting system, we must first find the 
interacting vacuum state |$). Our approach will consist in determining a unitary transfor- 
mation that maps the basic operators Ck«, sz,k into a new set Ck«, s/,k such that |$) can be 
written as 

|*>= E C ai ,..., aN fl 1 ...p NaN \^), (24) 

«i=T.i 

where e pa \^) = hq a \ty) = 0, and the new local Fermi operators fi )Ct are related to §i t i by Eq. 
(0), (we denote by (x) a former equation (x) with all operators tilded). 

The transformation is chosen to be unitary to ensure that the new operators continue 
satisfying the basic, defining algebra (g)-(g). The relationship between the initial and trans- 
formed operators is, thus, 



L = e f cl a e- f , s l>k = e f ~s^e- f , f t = -f . (25) 



c 



Since we are interested, at this point, in determining a ground state that does not break 
any symmetry, we will demand that T preserve the charge, lattice translational, and spin 
rotational symmetries. The simplest operator involving fermions that satisfies this condition 
is 

T = /2 f(k, k') [(c k |Ck'T - Ck^k'jSo.k'-k + V^CkTCk'^-Lk'-k + V^cl^y^s^-u , (26) 
where 

f *(k', k) = -f (k, k') (27) 

in order to satisfy = — T. Apart from this condition, the function T(k, k') is arbitrary, 
and it will be determined by demanding that the state (^) be the actual interacting vacuum 
state. 



For small values of J, the new operators are expected to be slight deformations of the 
initial ones. Thus, we will expand Eq. fl25|) in powers of T, retaining the first order terms 
only: 



J _ 7,\ 



°L + [ T > 6 Ll > s iM = + [T, s Jik ] . 



(28) 



These equations are the basis of the forecoming calculations. Using Eqs. (§) and (|l 
we can write them in a more explicit form: 

J 



«(kH4m + 7^(k',k) 



+ ( — ) c k'T(|) s °' k - k ' + ^ // 2^k'i(T)^ 1 ( _1 )' k-k ' 



^ k TU) ^ 4i\r 1 /2 

''km . - : f> (k) <{ PkTU) + 7^T7^*( k '' k ) [+(-)ck'T(i)So,k'-k + V^^i(T)S_i( 1 ),k'-k 



So,l 



,k = s ,k + — T(ki,k 2 ) [V24 lT c ka iS-i,k a -ki+k - V / 24 1 iPkaT*i,k a -ki+k 



J ~ 



si(-i),k = si(_i), k + (-)— T(ki, k 2 ) (c klT c k2T - c kii c k2i )si(_i) jk2 _ kl+k 



st 



C kiT(l) Ck 2i(T) S 0,k 2 -ki+k 



(29) 
(30) 
(31) 

(32) 



where the coefficient <S(k) in Eqs. fl29|) and ( |30"D is, for the moment, equal to one. 

With these equations, we can readily express the KLH ([16]) in terms of the transformed 
operators e pa , h m , and Si ti . To carry out this program, we have to: 



in the expression of the Hamiltonian 



i. Substitute Eqs. 



ii. write the c kQ , operators in terms of e pa and h m using Eq. ([I 



iii. normal-order these fermionic operators using the canonical anticommutation relations 

iv. and contract the su operators corresponding to the same lattice site using the multi- 



plication table (H) as in Eqs. (11) and (|T 



The execution of this program is straightforward but somewhat lengthy, since many terms 
are generated. To easily identify the procedence of each term in the following discussion, we 
will symbolically represent Eqs. (|29|)-(p2[) by = & + 5$, c = c + 8c, and s = s + 5s. 
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The kinetic and Kondo parts of H have the structures c^c and c^cs respectively. We 
will say, for example, that a particular term comes from (c^|<5c|s), meaning that it has been 
generated by executing the steps (ii)-(iv) on the part of the development of c'cs with that 
structure. 

We will now determine T(k, k'). Writing T in Eq. with the electron and hole 
operators (|17|) , we get that T = T\ + T 2 , where T\ and T 2 have the structures 



Ti = TTTTn q)e^-qV P + T(q, p)^- q e p Sp_ q l , (33) 



T, 



AN 1 / 2 
J 



2 



T (P, P') e J> W-P ~ T (q> qO^-q'^-qV-ql • (34) 



4ATV2 

The unitary transformation e T can be written, to leading order in J, as the product of 
e T2 and e Tl . Since T^l^) = 0, e T2 does not contribute to change the vacuum state |$), and 
it can be taken to be the identity without loss of generality: 

T(p,p')=T(q,q') = Vp, p', q, q' . (35) 

The presence of terms of the form eWs in H would imply that |$) is not an eigenstate 
and, hence, not the ground state of the system. In fact, by demanding the cancellation of 
these terms, we will be able to determine T(p, q). To leading order in J, this condition is 
given by 



(22) - T(p, -q) (e p -e q )x© = 0, (36) 



where the first term comes from (c^|c|s) and the second from (c^|5c) and (5c^|c). This 
condition and Eq. ( ^7|) imply that 

f (p, q) = -f (q, p) = . (37) 

Since e p > and £ q < 0, the denominator in this equation is always positive. However, 
when p and q are both very close to the Fermi surface (low-energy region), this denominator 
tends to vanish and T(p, q) becomes very large. This means that, in this region, higher order 
terms in J can give important contributions to e^h^s. A careful inspection of these terms 
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leads to the conclusion that their dominant contribution at low energies comes from (5ct|5c|s) 
and is given by 



3J^ 
16 



knp'q) + k%q') 

p' q' 



f(p,q) x (22). 



(38) 



Assuming a constant density of states (D = 1), the sums in this equation can be straight- 
forwardly calculated: 



-ET(p / ,q)= /„ ^-^~ln — , (-e q «l/2) 



N 



£p' - £q 



2e c 



(39) 



p' - ~v - H — Q 

Thus, the terms (|38|) diverge as e p and e q vanish and, consequently, the cancellation of 
the contributions to e)h)s is not achieved. 

This problematic state-of-affairs can be overcome by introducing a low-energy regulator 
in Eq. (|3~7|), namely, 

1 



T(p,q) 



-T(q,P) 



(40) 



where fj is a small positive function. Since at high energies this function is neglible compared 
to e p and — £ q , it will only play a role at very low energies and, therefore, we can approxi- 
mately take it as a constant equal to its zero-energy value. Obviously, fj will be determined 
by demanding that the sum of all contributions to the terms eWs [(pB]) and (|38[)] vanish 
when e p = e q = 0, namely, 



1 



3J\ 1 
In 



3?7 2fj 

Thus, fj is given, as a function of J, by 



x (22) 



0. 



(41) 



3J\ 1 
77 = In — 

' 8 2fj 



(42) 



Substituting (|40|), we see that the terms (^) are neglegible compared to the first order 
terms (|36| ) in the high-energy region. Since in this region ( |40| ) essentially coincides with 
(P?|), Eq. fl3"6D will continue to hold to leading order in J at high energies. In summary, with 
T(p, q) and fj given, respectively, by (140|) and (|42|), we accomplish our objective of cancelling 
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the terms eWs for all energies and, since if is a hermitian operator, this also implies the 
cancellation of the terms hes. 

To verify that unitarity is indeed preserved by this transformation, we can calculate the 
canonical properties (0) and (||) in terms of the transformed operators [using Eqs. fl29l)-(|32|). 
(|) and (§)], and check whether they continue to be satisfied. This calculation yields, for 
the terms proportional to 1, 



—2c 

•So,kSo,k' — ° Ok,-k' , 



(43) 
(44) 



where 
d(k) 



1 + — VT 2 (k',k) 

16iV y 



-1/2 



1 + 



3_J2 
16 



a 



N 2 



1-1/2 



£k| + V 

-1/2 



£ k 



1/2, 



-1/2 



1 + J 2 In — 

4:7] 



(45) 
(46) 



Since both <S(k) and a tend to 1 for small values of J, the canonical algebraic properties 
are essentially preserved. This would certainly not be the case if fj = and, thus, the 
introduction of a regulator is not only necessary to cancel the eWs terms at low energies, 
but also to achieve unitarity. The departure of d from 1 is very small and can be neglected. 
In the case of a(k), however, this departure is not that small at low energies. If kp is a wave 
vector of the Fermi surface, <5(kp) ~ 1 — 1/(4 In ^). For small J, the second term of this 
equation is certainly neglegible compared to 1, but it is large compared to J. Thus, to satisfy 
more accurately the canonical anticommutation relations for the two sets of operators, we 
will substitute the factor a(k) in Eqs. (^) and (|30D (which up till now was set equal to one) 
by the function ((45|). 

Having determined completely the relationship between the initial and transformed op- 
erators [Eqs. (p!^)-(|32|), (|40D, (|42|), and (f45|)1, we can now proceed to calculate H in terms 
of e pa , h qa , and §i ti following the steps (i)-(iv) detailed after Eq. (|32|). From the many 
contributions to each particular term, we shall only retain the dominant ones for small J. 
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Next, we summarize the results obtained, indicating where the dominant contributions stem 
from: 

• Hkmetic- The main contributions to the terms e^e and h)h come from (c^c), (6$\6c), 
(c^\8c\s), and (5c^|c|s), and are given by 

Hkinetic = E ( £ p)^pa^pa + E & '(~ 6 'q)^qoAq« > ( 47 ) 



p,« q,a 



where 



3 J 2 

%p) = " ^ E [(^P« 2 (P) - ^« 2 (q)) T(p, q) - 2a(p)a(q)] T(p, q) . (48) 



Assuming a constant density of states (D = 1), this equation can be developed into 



%) = e + 



3J^ 
16 



3J 2 1 \ £ + 1/2 7] 



1 - -7, In '— + 



16 e + fj J e + fj e + rj 



(s>0). (49) 



C: The essential contribution to the terms proportional to 1 comes from normal- 
ordering the operator — J2q,a H(— £ q)^qo^qa) which yields 

3 J 2 

q 8 

where -E$ = I] q 2e q = — iV/4 is the energy of the non-interacting vacuum state |$ )- 

HKondo'- The terms eWs anc l its hermitian conjugates foes are obviously absent from 
H, since this was the condition used to determine T(p, q). For the terms e^es and 
h'hs, the main contributions come from (c^|c|s) and are given by 

HKondo = a(p)a(p') x (20) + d(q)d(q') x (0) . (51) 

Hrkky'- The dominant contributions to the terms s^sj/j (i ^ j) are of order J 2 and 
come from (8&\8c), (&\8c\s), and (<5c^|c|s). This RKKY interaction is given by 

Hrkky = ^ E Jrkky(Rh — Ry)SjSj , (52) 
14 



where Sj are the vector operators related to by Eq. (0), and Jrkky(R) is an 
inter-site exchange coupling given by 

Jrkky(R) = ^ E [e P 3 2 (p)T(p, q) - 2a(q)fi(p)] f (p, q) cos[(q - p)R] . (53) 
p,q 

or, with the approximation a(k) ~ 1, 

W(H) = - 2) "K" - P» R 1 , (54) 

where F is the volume of the first Brillouin zone. It is clear from this equation that, 
in order to determine Jrkky(R), it is not enough to know the density of states, but 
we have to specify the entire function e^. In one dimension, however, W = D = 1 
and Ek = £-k imply that Ek = a\k\/n — 1/2 where a is the real lattice constant. This 



dispersion relation leads to 

J2 -1/2 -1/2 / 

Jrkky{R) = Y Jo dx J dy 1 



y 2 x 



1/2 + y — x + 1] 



cos[(l/2 + y- x)nR/a] + cos[(l/2 + y + x)irR/a] 
X I/2 + 2/-X + 77 ' ^ ^ 

where, since the integrals converge as fj — > 0, we can make the approximation of 
taking 77 = 0. From Eq. (0), it is clear that JijKR-y(R) is a long-range, oscillating 
coupling. Carrying out the integrals ( [55]) we find that, for the one-dimensional case, 
Jrkky{R — ifl) — 0.26J 2 , which shows that the exchange coupling between nearest 
neighbours is antiferromagnetic. 

(5$\5c), ($\6c\s), and (5c^|c|s) also produce terms of the form c'cSiSj (i 7^ j). The 
e'WsiSj parts and its hermitian conjugates should be cancelled by the J 2 terms of T 
in the same way as the leading orders of this transformation cancel the eWs and hes 
parts. The remaining terms, e'e§iSj and h^hSiSj, are associated with the interaction 
between charge carriers and local moments, and it may seem at first sight that they 
can be neglected because this interaction is dominated by HkoucLo- Nevertheless, in the 
presence of spin correlations without magnetic order (spin liquid), these terms would 
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give leading contributions to the low-energy electronic dispersion relation because the 
main contribution of the Kondo interaction vanishes ((<&|sj i|<&) = 0). However, since 
the effect of the spin correlations on the charged excitations will not be considered in 
this article, we will not include these terms in H. 

By carefully inspecting the rest of the terms generated by this change of coordinates, it 
can be shown that their effect is negligible compared to the interactions listed above. Thus, 
the expression of the KLH in terms of the new canonical operators e, h, and s is essentially 
given by 

H = C + Hkinetic + Hxondo + HrKKY ■ (56) 



III. SPIN DYNAMICS, MAGNETIC ORDER, AND ELEMENTARY 

EXCITATIONS 

On the basis of Eq. (|56]) we can now address the problem of magnetic order. To com- 
pletely specify the vacuum state |<E>), we still have to determine the coefficients C aij ... jCW of 
Eq. fl24]) which describe the correlations among the local / spins. Since all the terms of H 
containing tilded fermionic operators annihilate |$), these correlations arise exclusively as a 
consequence of the induced RKKY interaction. Thus, to determine the coefficients C ai) ... )QW 
we must demand the minimization of 

($\HrKKy\$) = ^JRKKY(B*-Bi)($\SiS j \*)- ( 57 ) 

#i 

This quantity [which would be zero if C ai) ... )CW were random coefficients] is negative, and it 
represents the energy decrease due to the formation of spin correlations. The total ground 
state energy is given by 

372 

E$ = E* -— Nln2 + (<S>\H RKKY \$), (58) 
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where E§ is the energy of the non-interacting Fermi sea, and the second and third terms 
represent the energy decreases due to the direct Kondo interaction and the induced spin 
correlations respectively. 

Thus, within the approach presented in this paper, the determination of the magnetic 
structure of the vacuum state for small couplings is reduced to finding the ground state of 
a Heisenberg Hamiltonian with an RKKY coupling given by Eq. fl54]) . This is, of course, a 
very hard task which we do not intend to undertake in this article. However, it should be 
noted that the situation corresponding to Hrkky is fundamentally different from that of 
a typical Heisenberg Hamiltonian coupling only nearest neighbours, because a long-range, 
oscillating coupling like ([54]) will inevitably lead to strong frustration and spin fluctuations, 
which tend to inhibit magnetic order and to favor the formation of a singlet RVB-like spin 
liquidjllil! or an antiferromagnet with strongly quenched local moments. 

The magnetic nature of the ground state has been usually attributed to the competi- 
tion between the Kondo and RKKY interactions, establishing a close parallelism with the 
discussion of the two-impurity case:i For weak couplings, the RKKY interaction (which is 
proportional to J 2 ) dominates and the two impurities order magnetically, while for stronger 
couplings, the formation of a Kondo singlet at each impurity site [with binding energy pro- 
portional to exp(— 1/ DJ)] is energetically favored and magnetic order is destroyed. Between 
these two extreme cases, we would have intermediate situations in which the local moments 
are only partially compensated. 

As we mentioned in the introduction, the extension of this approach from the two- 
impurity case to the Kondo lattice is very unclear due to the Nozieres exhaustion problem. 
Another objection that can be raised is that this mechanism of forming Kondo singlets to 
compensate the local moments would not be valid for J < and, therefore, could not explain 
why systems with ferromagnetic couplings like CeAl 3 do not exhibit appreciable magnetic 
order. 

The first and most popular attempt to implement in the Kondo lattice this idea of ex- 
plaining magnetic order as the result of the competition between the RKKY interaction and 
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the Kondo singlet formation was carried out by Doniach,El who studied a one-dimensional 
analog of the KLM (the so-called 'Kondo necklace') in which the conduction electrons are 
represented by pseudo spins. He took as variational antiferromagnetic trial function, the 
direct product of local wave functions consisting in a linear combination of the singlet and 
triplet states that can be formed with the local spins and pseudo spins. As a result, he found 
that magnetic order was destroyed for J > W. Although, as he argued, this should only be 
understood as a qualitative result, it is quite revealing since it essentially says that, if the 
formation of Kondo singlets was the mechanism responsible for the destruction of magnetic 
order, this would take place when J is strong enough to overcome the exhaustion problem, 
enabling the participation of the entire conduction band in the formation of Kondo singlets. 

Another well-known approach to the problem of magnetic order in HF systems is due 
to Coleman and Andrei0, who claimed that the Kondo scattering of low-energy conduction 
electrons and local moments could stabilize the spin liquid against an antiferromagnetic 
state which is supposed to have a rather similar energy. This idea is still being investigated 
by some authors.il 

In contrast to these approaches, we have found in the last section that when the RKKY 
interaction is explicitly extracted from the Hamiltonian, the Kondo interaction is modified 
in such a way that it no longer influences the ground state. Thus, as we mentioned, the 
formation of a spin liquid has to be attributed to the long-range character of the RKKY 
interaction which produces strong frustration and enhances spin fluctuations in the system. 
Of course, this picture holds only as long as the ground state |$) remains stable. In fact, we 
will later find that an electronic instability appears for stronger couplings and, at that point, 
the Kondo coupling between conduction electrons and local moments will indeed influence 
the spin correlations. 

It should be noted that, since Hrkky is invariant under the change of J by —J, the 
discussion of magnetism is (before the occurrence of the mentioned instability) the same for 
both ferro and antiferromagnetic couplings. Finally, we would like to point out that other 
effects not contemplated in the KLM like the existence of a small f-charge fluctuation could 
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also contribute to the stabilization of the spin liquid. 

The approach presented in this paper, also leads to a physical interpretation of the 
elementary excitations of the system and their role in the various observable quantities 
which is fundamentally different from the most popular interpretations. We have seen in 
the previous section that the expression of the KLH in terms of electron and hole operators 
associated with the actual vacuum state is inevitably accompanied by the appearance of 
an explicit RKKY interaction which leads to the formation of spin correlations. Thus, 
the physical picture of the system at zero temperature motivated by our formalism is that 
of a backgroung of spin correlations (a spin liquid) in which two kind of modes can be 
excited: neutral spin modes associated with deformations of the spin liquid and charged 
modes corresponding to the excitation of electrons and holes in the system. 

The accurate determination of the state of the spin system and its excitations would 
require finding the ground state and spectrum of Hrkky which, as we mentioned, is beyond 
the scope of this paper. Nevertheless, we can still draw an important conlusion without 
having to solve this Hamiltonian: The thermal breakdown of the spin correlations will 
produce an entropy increase of the order of iVln2; since H^kky is a very weak, second 
order interaction, the characteristic temperature T^l associated with the spin correlations 
will be very small and, therefore, this enormous entropy increase will take place in a very 
small temperature interval, leading to a huge specific heat of spin origin. 

Thus, as we mentioned in the introduction, the formalism developed in this paper sup- 
ports the thesis of Kagan, Kikoin, and Prokof 'ev,0 who argued that a strong mass renormal- 
ization of the low-energy conduction electrons cannot explain the high value and universal 
character of the entropy measured at Tkl, concluding that the heavy particles responsible 
for the low-temperature properties of the Kondo lattice should be neutral fermions of spin 
origin, as opposed to the charged heavy fermions of the slave-bosonlil'0 and variational^ 
approaches. 

In the next two sections, we study the charged excitations of the system. 
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IV. CHARGED EXCITATIONS: PERTURVATIVE MODES. 



Since we are interested in studying the zero-temperature charged modes for a Kondo 
lattice in the normal state, we shall assume that |$) is a singlet spin liquid, namely, 

S|$) = 0, (59) 

where S is the generator of spin rotations. The spin rotational symmetry will be very useful 
to characterize the states as eigenvectors of S 2 and S 2 with definite eigenvalues. For a general 
state of the form V|$), where V is a linear combination of products of e^ a , h^, and 
operators, we can write e = e iaS Ve~ 1 ' aS \§), which implies that the action of S on a 

state is given by 

sy|$) = [s,v]|$>, (60) 



We would like to emphasize that, since the transformation (26) preserves the spin symmetry, 
the commutation rules of S with the initial and transformed operators are exactly the same. 
In particular, Eq. (|60| ) implies that the states ej, a |<&) and s/,k|3>) form, respectively, an 
s = 1/2 charged doublet and an s = 1 triplet of spin fluctuations. The fact that S z §i j i c \$>) = 
^z,k|$) justifies the notation introduced in Eq. (|7|) for the s operators. 

Aside from the contributions of spin correlations, the dispersion relation corresponding 
to the states ej, a |$) and ^q Q |$) is given by Eq. ( fl9|) and contains a gap of magnitude 
2E(0) = fj + 3J 2 /16. This statement, however, does not have any physical significance 
because these states are not the charged excitations of the system. 

For small J, the actual charged modes are expected to be slight modifications of the states 
ej^J^) and /i qc J$) and, consequently, we will refer to them as perturbative modes. However, 
as we argued in the introduction in relation to the Nozieres exhaustion problem, we should 
also consider the possibility that the Kondo interaction would produce a strong correlation 
between low-energy electrons and spin waves, leading to the formation of collective states. 
These strongly correlated modes will be investigated in the next section. 
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When we determined e T in Section II, we argued that the transformation e T2 [Eq. 
did not change |<&) and, therefore, could be set equal to the identity. Nevertheless, since e T2 
does change the basic operators Ck a and s^k, we can use the degrees of freedom associated 
to this transformation to get a new set of canonical operators Ck a and s^k which still satisfy 
e p «|$) = /iqa|$) = and which are optimal in the sense that e pc J$) and ^q Q |$) are the 
best possible approximations to the actual charged modes. 

In the next section, we will calculate this second transformation (T), and will find that 
Cka, si,k are slight deformations of the tilded operators. Since Ck a , s^k are optimal operators, 
they should be the basis to construct the actual charged modes. In fact, they will be 
used in the calculation of the strongly correlated modes. To construct the perturbative 
modes, however, we can also use the tilded operators because the composition of two small 
deformations is still a small deformation. In fact, the direct calculation with the tilded 
operators is more reliable because we avoid the inevitable approximations associated with 
carrying out a second transformation. 

Since no symmetry is broken and no bound states are created at this point, there is no 
physical reason for the appearance of a gap in the charged spectrum. In fact, we will now 
calculate the dispersion relation for the perturbative modes to leading order in J neglecting 
the effect of the spin correlations, showing how the gap is indeed closed in this order. The 
perturbative modes must have s — 1/2. Due to the spin rotational symmetry, the two states 
in a spin doublet will have the same energy and, therefore, we can select, for instance, the 
state with s z = 1/2 to carry out our calculations. Moreover, since the situation considered 
is symmetric with respect to electrons and holes, we only need to calculate the spectrum 
corresponding, for example, to the electrons. It can be seen, using Eq. (0), that a small 
deformation of the state e p J<I>) having negative charge, s z = s = 1/2, and wave vector p, 



has the structure 

|P£ pT > = a(p) [e pT - i ^/3 p (p , )(4'T S o lP -p' + V^'i^p-P')] |$> , (61) 
where PE stands for 'perturbative electron', /3 p (p') is an arbitrary variational function, and 
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a(p) is a normalization factor given by 

3J 2 



a{p) 



1 + 



i6r 



-1/2 



(62) 



The function /? p (p') should be selected to turn the state ( |61| ) into an energy eigenstate. 
This condition can be implemented by demanding that /3 p (p') be an extreme of the functional 

Obviously, to calculate this quantity, we will use the expression of H in terms of the tilded 
operators [Eq. (|56|) 1. The scalar product in this equation can be developed with Eqs. (§), 
(|), and (|) into an expression in which only scalar products of the form ($1^,^ . . . Si m ,i m \®) 
with ii, . . . ,i m corresponding to different lattice sites, leave to be determined. These are 
precisely the spin correlation functions, and depend on the magnetic structure of |$) arising 
from the RKKY interaction. As we mentioned, the calculation of these correlations is beyond 
the scope of this paper and, therefore, its effect on the charged modes of the system will 
be neglected. Thus, we will calculate these modes assuming a disordered spin lattice, which 
means that all the spin correlation functions vanish, enabling a simple computation of the 
scalar products. Since the constant C in Eq. (|56| ) coincides, in this case, with the ground 
state energy, we will remove it from H assigning, in this way, zero energy to |$). This 
means that the energy of every state is measured with respect to the vacuum state, and 
represents, therefore, its excitation energy. Obviously, this energy has to be positive; in 
fact, the existence of a state with negative energy, would imply an instabily of |$) and a 
transition to a new phase. 

With all these conditions taken into account, the calculation of E([j3],p) yields 

3J 2 



E 



E e T 



EM 



16 

-E(e 



a 2 (p) I a[£p) 



r 



dp'a(e p ,) /3 P (p')+/W) 



(63) 



J 



/3p(p')A>(p') + TTp / rf Pl rf P2«(£pi)«(£p 2 )/3p(Pl)A>(P2' 



where a(e) and E(e) are given by Eqs. (|45D and (|49| ) respectively. The last term of this 
equation is of higher order in J than the others, and can be neglected. 
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If the gap is indeed closed, E(0) = E(\p],pp) (where pf is a wave vector of the Fermi 
surface) must vanish. This quantity can be written as E(0) = E(0) — A, where A represents 
the energy decrease due to the Kondo interaction, and is given by 

A = 3JMp F ) J dp | 5(0)fi(£p) r (p) + (p)] _ ^(s p ) - E(0)} /?*(p)/?(p)} , (64) 

where /3(p) = /3 PF (p). Functionally differentiating with respect to /3(p) in this equation, we 
find that the condition of extreme of E(0) is given by 

£7(e p ) - -E(O) + A 

and substituting this expression in Eq. (0), we obtain the following equation for A: 

A = a 2 (0) / de^ . (66) 



16 w Jo E(e) - E(0) + A 
As we already mentioned, a(e) can be substituted, to leading order in J, by 1. Also, since 
E(e) — E(0) vanishes at e — 0, we can approximate it by e. Thus, with the approximations 

5(e) ~ 1 , E(e) ~ E(0)+e, (67) 

Eq. ( |66D is turned into 

A = ^l n J_ (68) 

16 2A' v ; 

which implies that A = fj/2 + 0(J 2 ). Since E(0) = r//2 + 3J 2 /32, this proves the cancellation 
of the leading contributions to E(0). We will not verify the cancellation of the J 2 terms 
because, since these are lower order contributions, to accurately calculate them, we should 
have to control all the approximations that have enabled us to carry out our computations 
analytically. Since, as we already mentioned, there is no physical reason for the formation 
of a gap, we will assume that the J 2 contributions cancel as well. 

Finally, we would like to point out that if the spin correlations were considered, we should 
include in H the terms e^e§i§j and h^hSiSj {i ^ j) that we mentioned at the end of Section 
II because, in this case, they give leading contributions to E(0) which compete with those 
steming from the kinetic and Kondo terms. 
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V. SECOND TRANSFORMATION, STRONGLY CORRELATED CHARGED 
MODES, AND VACUUM INSTABILITY 



In this section, we shall investigate the possible formation of strongly correlated charged 
modes. As we already argued in Section IV (and will become even more clear at the end of 
this section), these modes should be written in terms of optimal fermionic operators repre- 
senting the best possible approximations to the actual (perturbative) charged excitations. 
Thus, we will start by carrying out a second unitary transformation to determine these op- 
erators, and will find the expression of the Hamiltonian in terms of them. Then, with these 
building blocks, we will construct the strongly correlated charged modes and will discuss 
their nature in the case of ferro and antiferro magnetic couplings. Finally, we will show that, 
for sufficiently strong couplings, these modes condensate in the vacuum state leading to an 
instability. 

In analogy with the notation (x), (x) will denote the equation equivalent to (x) for the 
optimal operators Ck«, sz,k- Since we want to continue preserving all the symmetries of H, 
the second unitary transformation T that we will now determine must have the structure 
([26|), and the relations between the tilded and hatted operators will be given by Eqs. fl29|)- 
(|32|) with d(k) equal, for the moment, to one. To ensure that e pQ |<3>) = /iq a |$) = 0, e T must 
leave |<3>) invariant, which requires that T(p, q) = T(q, p) = 0. 

We shall take as starting point in the following discussion, the expression (|56|) of H in 
terms of the tilded operators with the approximations (^)7|). The use of these approximations 
in the calculation of quantities in which the dominant contributions cancel will be justified. 
The main effect of the Kondo interaction is, obviously, to renormalize the charged excitations. 
Since e^J^) and /i q J$) are optimal approximations to the perturbative charged modes, 
they will include these renormalization effects and, therefore, they will be subject only to a 
small residual Kondo interaction. Thus, when we write H in terms of the hatted operators, 



the bulk of the Kondo interaction should disappear. In analogy with equation (|36D, the 
condition for the cancellation of the terms e^es and Whs to leading order in J is given by 
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T(p,p') 



r(q, q') 



(69) 



Again, it can be shown that the singularities of these expressions at the points where 
£k = £k' lead to a loss of unitarity. Thus, we need, once more, to introduce a regulator to 
tame the singularities. In contast to e p — £ q [Eq. which is always positive, e p — e p i and 
£q — ^q' can have positive and negative values and, therefore, to avoid the vanishing of the 
denominators of Eq. (|69|), we should add to them a small imaginary quantity. Thus, the 
regulated forms of T(p, p') and T(q, q') can be written as 

1 1 



T(p, p') 



— + — 

Bp - £ P ' +iV Bp- Bp' - trj 



Bp') 2 + fj 2 



(70) 



T(q, q') 



(71) 



(£q-£q0 2 +V 2 ' 

where, to re-establish unitarity for small couplings, the real constant fj should vanish more 
slowly than J 2 as J — > 0. Again, as we did with the first transformation in Section II, to 
attain unitarity to the next-to-leading order in J, we should introduce in Eqs. Q2§D and ( ^U[ ) 
the factor d(k) given by 



a(k) 



1 + ^yf 2 (k',k) 

16A^ . 



-1/2 


3 T 2 r 




1 + 




16 io 



dx- 



(kk| 



-1/2 



(72) 



kk| - x) 2 + 17 2 ] 

To completely specify this second transformation, we just need to determine fj. We saw in 
the last section that the Kondo interaction lowered the energy of the charged states around 
the Fermi level; in fact, the equations (|65| ) and (|66|) which determined the state \PE PF ^), 
are nothing but the conditions for minimizing the energy of this state. Since e|, f|<I>) must 
be the best possible approximation to \PE PF -^), we will determine fj by demanding that its 
energy be minimized. 

If the spin correlations are disregarded, the energy of the states ep a |$) is given by the 
coefficient E(p) of the kinetic term (J2 pa ^(p)^ Pa ^pa) °f H expressed in terms of the hatted 
operators. In particular, e J$) has the energy E(0) = E(0) — A where A is given, with 



the approximations 



for a(e) and E(e), by 
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3J 2 fl/2 , 2 , £ 2 



A = — / de h^OHs)-^)^— - ^— -, (73) 

where the function a(e) is given by Eq. (|T2|). After some calculations, it can be shown that 
the condition of extreme of A (dK/dfj = 0) leads to the following equation for r): 

77 = In — , (74) 

' 64 2))' K J 

and that the corresponding (maximum) value of A is given, to order J 2 , by 

A = ^( ln ±_iV (75) 
16 V 2f) 2 J K J 

E(0) can also be written as E(0) = E(0) +A — A = A — A. Since the dominant parts of A 
and A cancel in this equation, E(0) will be proportional to J 2 and it seems that, to calculate 
the coefficient of proportionality, an accurate computation controlling all the approximations 
[like (0)] would be needed. However, although these approximations modify the order J 2 of 
the energy decreases A and A corresponding to the states \PE PF ^) and ej, J$), their effect 
on these two states must be very similar and, therefore, they should not essentially change 
E(0) = A — A. Thus, by substituting in this equation the values of A and A calculated with 
the same approximations in Eqs. ( j68| ) and (|75|), we should get an estimate of E(0). The 
result obtained in this way is 

3 T 2 / 7T 1 \ 

E(0) = A-A = — (ln- + -) = 0.048J 2 . (76) 



16 V 4 2, 

Having determined fj [Eq. ([74])], the relations ([2^)- (|32|) between the tilded and hatted 



operators are completely specified. In analogy with what we did in Section II with the 
first transformation, by substituting these relations in Eq. fl56|), normal-ordering the hatted 
operators corresponding to electrons and holes, and contracting the operators s^j at the same 
lattice sites, we obtain the expression of H in terms of the optimal coordinates e pa , hq a , and 
s«,k- By carefully examining all the terms generated in this process, we can discriminate the 
ones which are dominant for small J. After this analysis, we find that the new expression 
of H is essentially given by 
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H — C + H kinetic + H Kondo + H RKKY ■ (77) 

where 

Hkinetic = E E(e p )el a e pa + £ £(-e q )fct A qa with %) ~ £(0) + £ , (78) 



= (e „ _ £ S + y x (2 " 0) + fa - v)' + y " ^ ' (79) 

Hrkky = - E ■^flKK'y(R^ — li/ )S,S ; . (80) 

Thus, as expected, the main effects of the second transformation are to strongly reduce 
the gap in the kinetic term and to substitute the Kondo coupling by a residual Kondo 
interaction which is effective only between electrons or holes having very similar energies. 
In contrast to what happened with the first transformation, the technical reason for the 
appearance of this residual Kondo interaction is that, in this case, the contributions of 
higher order terms are not important and cannot cancel the Kondo coupling in the regions 
where ~ 

Armed with this (optimal) expression of H, we will now look for strongly correlated 
charged states, namely, we will investigate the possibility that an electron or a hole combines 
with a wave of spin fluctuations forming a collective state. Again, since the situation is the 
same for positively and negatively charged modes, we will only need to investigate, for 
example, the last ones. The structure of such a correlated state having negative charge, 
wave vector k and s = s z = 1/2 must be 

\SCE^ 1/2 ) = N- l ' 2 AB( V ) (e^k-p + v/2e^s ljk _ p ) |$> , (81) 

where SCE stands for 'strongly correlated electron', B(p) is a variational function, and A 
is a normalization factor given by 



A 



d P B*(p)B(p) 



-1/2 



(82) 



Neglecting the spin correlations, the energy of this state, measured with respect to the 
ground state energy, is given by 
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Esc,i/2[B] = A 2 



^ JdpE(e p )B*(p)B(p)-^ Jdpdp'B*(p')J(e p ,,e p )B(p) 



(83) 



where J{e\,E2) = Jrj 2 /[(ei — E2) 2 + f} 2 } is the residual Kondo coupling of Eq. (|79|). 

Obviously, the variational function B(p) will be determined by demanding that E,sc,i/2 
is minimized. This condition is given by 

{E{e p ) - E sc , 1/2 )B(p) - J dp'J(E p , E p/ )B(p') = . (84) 

From Eqs. ([83] ) and (|84]) , we should be able to calculate Esc,i/2 and B(p). It should be 
noted, however, that, to minimize Esc,i/2, the kinetic part in Eq. ([33] ) must be kept small, 
which requires that B(p) have all its weight concentrated in a very narrow layer around the 
Fermi level. Since the residual Kondo coupling is essentially constant and equal to J in this 
region, we will make the approximation of substituting in Eq. (]8~4"D the function J(£ p ,£ p r) 
by J, and will prove later on that this is indeed a valid approximation. Thus, from Eq. (|84]) 
we have 

B(p) = g- - Z p , Z = ^ / rf P 5(p) • (85) 

Substituting the first of these equations into the second one, we get the following condition 
for E S c,i/2- 

1 = - ^ . (86) 

2 Jo e + E{0) -E sc ,i/2 

This equation has a solution only for antiferromagnetic couplings (J > 0), and it is given by 

e~ 2/J 

Esc,i/2 = E(0) — O1/2 , ^1/2 = — - — • (87) 

For ferromagnetic couplings, the local moments and the spin of the conduction electrons 
will tend to align parallel to each other and, if a collective state is formed, it should have 
the structure 

\SCE K3/2 ) = N-^CDip^^-pl®) , (88) 
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which is the general form of a state composed of a charge carrier and a spin wave having 
wave vector k, negative charge, and s = s z = 3/2. Again, D(p) is a variational function 
and C a normalization factor. The energy of this state in the absence of spin correlations is 
given by 

If-, 1 



E SC ,3/2[D]=C 2 [ f JdpE(e p )D*(p)D(p) + — Jdpdp'D*{p') J( £p ,, e p )D(p) 



(89) 



Following the same steps used to estimate E SC) i/2, we arrive to an equation for E S c,3/2, 

i=4 r , de , (90) 

4 Jo £ + E(0)-Esc,3/2 
which has solution only for ferromagnetic couplings and it is given by 

e 4/J 

EsC,3/2 = E(0) - Q 3 / 2 , ^3/2 = —y- ■ (91) 

Thus, summarizing, we have found that in ferro and antiferromagnetic Kondo lattices 
there is a tendency to form collective states of spin 3/2 and 1/2 respectively, and that, for 
couplings with the same strength, the formation of these states is much more favored in the 
antiferromagnetic case (f^i/2 ^> ^3/2)- 

If fl stands for either Qi/ 2 or Q3/2, the probability of finding in a SC mode an electron with 
wave vector p in the layer e < e p < e + de is P(e)de, where P{e) ~ Q/(e + Q) 2 . This means, 
for instance, that 50% of the electrons in these modes have their wave vectors in the layer 
< e p < £1. Thus, if f2 <C E(0), the SC modes cannot form because it is energetically more 
favorable for the electrons in this layer to remain in the PE modes. However, if Q is smaller 
but of the order of -E'(O), the low-energy electrons are energetically allowed to participate 
in the formation of both PE and SCE and this will result in a very strong renormalization 
of the electron (and hole) masses around the Fermi level. Finally, if Q > E(0), the energy 
of the SCE will become negative, which means that these states condensate in the ground 
state producing an instability and the transition to a new phase. 

It should be remembered that, to simplify the calculations, we took a constant density 
of states D — 1. However, in a real band D(e) vanishes at the edges and it is enhanced in 
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the middle of the band (Fermi level). Since, as we have already mentioned, the SCE are 
composed of very soft electrons, only the density of states at the Fermi level (Dp) plays a 
role in their formation. Thus, more realistic expressions of Q are 

e -2/JD F e 4/JD F 
= , ^3/2 = ■ (92) 

Since the value of Q depends very strongly on the exponents of these expressions, it is very 
important to use the correct value of Dp when we try to estimate J c , namely, the value of 
J at the phase transition point (Q = E(0)). For example, in the antiferromagnetic case we 
have (in units of W) J C {D F = 1.3) = 0.34, J C (D F = 1.5) = 0.27, and J C (D F = 2) = 0.17. 

At this point, we would like to show, as we promised, that the approximation J(e p , e p /) ~ 
J used to determine E S c is indeed justified: The average value of the residual Kondo coupling 
j(£i,e 2 ) in a SCE state is given by 

(J) = f 1/2 [ 12 de l de 2 P(e 2 )J(e ll e 2 )P(e 1 ) . (93) 
Jo Jo 

Substituting the expressions for P(e) and J(£i,e 2 )) this equation can be developed into 

(J) „ fm dx fx + 2 2, , A ,„ A . 

J Jo x 2 + r\ l I \¥ \x + 1 x J 

Since the second term of this equation is small when fj/Q ^> 1 (the actual physical case), 
(J) ~ J, and our approximation is justified. For example, in the antiferromagnetic case 
with Dp = 1.5 we have, at J = J c , rj/Q = 8.6 and (J) = 0.81J. 

It is interesting to note that in our calculation of the SCE starting from the Hamiltonian 
(j56|), we have reduced the minimum energy E(0) of the states ej, w |<&) in two steps: First 
finding optimal operators e pa which strongly reduced the minimum energy from E(0) to 
E(0) using the high-energy tail of the Kondo coupling and, second, constructing with these 
operators SC modes which further reduced the energy using the low-energy part of the 
Kondo coupling. This makes it very clear why we needed to find the optimal operators to 
construct the SCE: If we would have built these modes with the operators e pa and s^k, we 
would not have cacelled the high-energy tail of the Kondo coupling, and the resulting states 
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could not be considered as proper states because they would continue interacting strongly 
with the lattice of local moments. 

As we mentioned, we have always neglected the influence of the spin correlations on 
the charged modes. This influence is actually expected to be important in the quantitative 
description of the SC modes. At first sight, it seems that these correlations would tend to 
inhibit the formations of SCE because an additional energy (proportional to J 2 ) would be 
needed to excite the spin waves present in these modes; however, the changes in the scalar 
product introduced by the spin correlations would also modify the contribution of the Kondo 
interaction to the energy and the overall effect is hard to predict. 

The most important challenge posed by this article is the description and characterization 
of the new phase. Obviously, since we know experimentally that the instabilities in HF 
systems are associated with antiferromagnetim and superconductivity, it is important to 
elucidate whether or not the instability that we have described is associated with the onset 
of superconductivity. In any case, it is clear that the condensation in the ground state of 
collective charged states will dramatically change the diamagnetic and transport properties 
of the system. It should also be noted that, since in these condensated states, spin waves 
strongly couple to conduction electrons, the transition to the new phase will be accompanied 
by changes in the spin correlations. Thus, since we have attributed the high values of the low- 
temperature specific heat to the thermal breakdown of soft spin correlations, an important 
change in the specific heat is expected to occur at the transition point. In fact, since the 
state of the spin system affects the SC modes and the condensation of these modes alters the 
spin system, an accurate calculation of the transition point would require a self-consistent 
calculation involving the two systems. 

VI. SUMMARY AND CONCLUSIONS 

In this article, we have assumed that the KLM incorporates the essential physics of HF 
systems and, investigating this model at zero temperature and half filling with a 'variational 
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perturbation theory', we have obtained the following physical picture: 

• The non-interacting system consists in an electronic Fermi sea with electrons and holes 
as elementary excitations, plus a completely degenerate (random) spin lattice. 

• For small couplings, the charged spectrum remains basically unchanged but an induced 
RKKY interaction breaks the degeneracy of the spin lattice. 

• At this stage, the RKKY interaction is the only one responsible for the spin dynam- 
ics. This very long-range, oscillating interaction implies strong frustration and the 
enhancement of spin fluctuations, which should explain the formation of a spin liquid 
or the appearance of a weak magnetic order depending on the crystal structure of the 
compound. 

• The very high values of the low-temperature specific heat measured in these systems 
are attributed to the enormous entropy increase associated to the thermal breakdown 
of soft, RKKY-induced spin correlations. 

• For higher values of the Kondo coupling, the low-energy electrons and holes are en- 
ergetically allowed to participate in the formation of both uncorrelated and strongly 
correlated charged modes which results in a very strong renormalization of the low- 
energy charged spectrum. The formation of strongly correlated charged modes is much 
more favored in antiferromagnetic Kondo lattices than in the ferromagnetic ones. 

• Beyond a critical coupling J c , the strongly correlated modes condensate in the ground 
state producing an instability and a transition to a new phase. This transition must 
be accompanied by an important change in the specific heat because the association of 
spin waves to conduction electrons to form collective states will modify the localized 
spin system. 

The main future challenges posed by this article are the detailed study of the spin 
correlations, its inclusion in the calculation of the charged modes, and the description and 
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characterization of the new phase. 
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